A Dynamic Approach to the Thermodynamics of Superdiffusion 
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We address the problem of relating thermodynamics to 
mechanics in the case of microscopic dynamics without a finite 
time scale. The solution is obtained by expressing the TsaUis 
entropic index q as a function of the Levy index a, and using 
dynamical rather than probabilistic arguments. 
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The problem of establishing a connection between dy- 
namics and thermodynamics is still the object of wide 
interest among physicists. Recently rather than deriving 
Brownian diffusion from statistical mechanics, and con- 
sequently from thermodynamics, the inverse procedure 
was applied, of derivingBrownian motion from merely 
dynamical arguments . It was also shown that this is 
equivalent to deriving thermodynamics from mechanics. 
However, this mechanical derivation of thermodynamics 
rests on the assumption that a finite microscopic time 
scale exists. In the last few years an increasing interest 
has been devoted to generalized forms of Brownian mo- 
tion that are incompatible with the existence of a finite 
microscopic time scale, notably the processes of superdif- 
fusion, and especially the Levy diffusion [^-1). Conse- 
quently, the dynamic derivation of these processes, at a 
first sight, seems to be incompatible with thermodynam- 
ics. 

It is well known Q| that these diffusion processes, in 
the one-dimensional case, are characterized by probabil- 
ity distributions p{x, t) whose Fourier transform in the 
symmetric case reads 



p{k,t) = exp{~b\k\"t), 



(1) 



where a is the Levy index ranging, in principle, in the 
interval < a < 2, and b denotes the diffusion intensity. 
Note that according to the analysis of , the processes 
corresponding to the more restricted range 



I< a < 2 



(2) 



are compatible with dynamic derivations, including those 
of deterministic and Hamiltonian nature The work 
of Ref. ||] illustrates a dynamic derivation resting on a 
stationary correlation function of the microscopic fluctu- 
ation, ^1^(0 0^ 1/^'' with < (3 < 1, namely with an 
infinite time scale, thereby making it impossible to de- 
rive thermodynamics as in Ref. |^ . Here we show that in 
spite of this limitation, these dynamics are not incompat- 
ible with thermodynamics, provided that a non extensive 
rather than an extensive entropy is adopted. 



Note that the diffusing system of Ref. [g| has a fi- 
nite velocity with only two possible values, W and —W, 
corresponding to two distinct states of the system, and 
the time duration of these two states, t, is character- 
ized by the probability distribution ip{t). The function 
ip{t), in turn, is proportional to the second time deriva- 
tive of ^w{t) and consequently oc l/t^"*"" with 
1 < a < 2. A jump of a given length \x\, with prob- 
ability density n(x), is not instantaneously made, but it 
takes place through a time interval proportional to the 
distance |a;| travelled by the particle. Consequently, the 
correlation function with an infinite lifetime "^wlt) re- 
sults in H(x) oc l/|a;| 

A crucial aspect of this dynamic approach is that the 
distribution is characterized by finite moments, in spite 
of the fact that the moments of the Levy diffusion, from 
that of the second-order on, are divergent. The appar- 
ent conflict between these dynamically derived diffusion 
processes and the traditional Levy processes is settled 
by noticing that the dynamic derivation implies the ex- 
istence of a p ropagation front, characterized by ballistic 
peaks [^O|,0, and that the probability distribution of x 
between these peaks takes the shape of a truncated Levy 
process. The shape of an ideal Levy process is reached, 
instead, in the time asymptotic limit, as a consequence 
of the fact that the intensity of these peaks tends to zero, 
although very slowly. 

It has to be pointed out that the historical founda- 
tion of the Levy processes is probabilistic in nature . 
We note that, according to the perspective estabhshed 
by Jaynes there is a close connection between the 
second principle of thermodynamics and a probabilistic 
approach to statistical mechanics. This means that the 
thermodynamic foundation of the Levy diffusion should 
be more natural than its dynamic foundation. On the 
contrary, as recently pointed out by a group of authors 
|l^-p^, the moment divergencies raise problems also to 
this thermodynamic foundation. These authors proved 
that if a proper use is made of the TsaUis entropy ||l^ 



s,[nix)] 



l-/_°lH(x)^dx 

(1-9) 



(3) 



the challenge posed by the infinite moments can be suc- 
cessfully addressed, since the constraint on the moments 
is set using H(a;)' as a weighting function thereby ensur- 
ing convergence even where the standard weight would 
produce divergencies. Using this kind of weight, the Levy 
process can be derived by means of the standard proce- 
dure of entropy maximization resting on the technique 



1 



of Lagrange multipliers. Later on in this paper we shall 
give more details on this approach. For the time being, 
let us limit ourselves to noticing that the maximization 
of the non-extensive entropy of Eq. ^ led the authors 
of Refs. to the conclusion that the entropic index 

q is related to a as follows: 



(4) 
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To realize the main purpose of establishing a connec- 
tion between dynamics and thermodynamics, we must 
make the result of Eq. compatible with the first steps 
successfully made by the authors of Refs. [^-^ at the 
dynamic level. These authors laid the foundation for 
the connection between the entropic index q and fractal 
dynamics. It is argued Q that it is convenient to gen- 
eralize the definition of the Kolmogorov entropy, which 
is a sort of change per unit of time of the Gibbs entropy, 
by replacing the Gibbs entropy with the Tsallis entropy. 
The generalized form of Kolmogorov entropy, called Kg, 
is easily shown to yield the generalization of the well 
known Pesin identity pi| to 



Kq — \q, 



(5) 



with \q playing the role of a generalized Lyapunov co- 
efficient. This important result is obtained by study- 
ing the distance A?/(t) between two trajectories depart- 
ing from two distinct initial conditions, at the distance 
A2/(0) the one from the other, and assigning to 5{t) = 



lim 



, Ay(t) 
|Ay(0)HO I 



the following form: 



5{t)^{l + {l-q)\t)- 



(6) 



It is evident that the traditional exponential sensitivity to 
the initial conditions is recovered from Eq. by setting 
g = 1. Note that we adopt the symbol y to denote micro- 
scopic trajectories so as to avoid confusion with the prob- 
abilistic and macroscopic level of description, denoted by 
the symbol x. 

In Refs. 1^,^ it has been discussed how to derive su- 
perdiffusion on the basis of dynamic properties, with 
a procedure, valid also in the Hamiltonian case, which 
however recovers the results of the continuous random 
walk method of 11 2^,|2^ when applied to the same 
dynamic systems. These dynamic systems are bimodal 
maps, with two laminar regions, derived from that of 
Geisel and Thomae Q. With respect to the laboratory 
reference frame, the particles are characterized by the ve- 
locity W or —W, according to whether they move in the 
right or left laminar region. The analytical form of the 
map describing the motion in the left laminar region is 
given by [pi: 



Vn 



Vn + ay„ 



(7) 



with a > 0, 2; > 1. This map extends from y = Oioy — d 
with the parameter d defining the size of the laminar re- 
gion. The right laminar region, with the same size d, is 



separated from the left by a chaotic layer and its analyti- 
cal form is the same as that of Eq. with y denoting in 
this case the distance from the right border of the map. 
The particles starting from an initial condition very close 
to y = are assumed to be driven by the continuous time 
approximation to Eq. (Q): 



dy_ 

dt 



(8) 



Thus it is possible to evaluate the distribution of waiting 
times in the interval (0,c?). The assumption that trajec- 
tories are injected with uniform probability yields ||22| 



a[l-|-d^"-^a(z- 1)<]" 



(9) 



The right laminar region yields for the corresponding 
waiting time distribution the same analytical form as 
that of Eq. (^). Note that the research work of Refs. 
||7|,p|p^ has established the connection between the in- 
dex z and the Levy index a to be given by 



(10) 



z-1 



Note also that the value z = 2 is the border with the 
region where the first moment of tpit) diverges, and rep- 
resents the limit of validity of the dynamic treatment . 
For values of z < 3/2, on the contrary, the second mo- 
ment of ip{t), as well as the first, becomes convergent. 
Thus z — 3/2 is the border between the Levy and the 
Gaussian region. 

To assess the validity of the entropic prediction we 
calculate 6(t) using the dynamics of Eq. and we prove 



S{t)^{l-aiz-l)y'o-hy 



(11) 



where yo is the initial condition of one of the two trajec- 
tories under study. As done in this initial condition 
is set so close to y = as to make it possible to adopt the 
continuous time representation of Eq. (j^ . Furthermore 
the assumption is made that |Aj/„/^| << 1. Notice 
that in the numerical calculations 1 1 1 1 it is usually set 
d = 1/2. Thus the divergence of S{t) of Eq.([l^) occurs at 
times much larger than the escape time from the laminar 
region. By comparing Eq. (|ll| ) with the general property 
of Eq. ^ , and using also Eq. (|l^, we obtain 

2 + a 



q = 



1 



a 



and 



Xq = a- 



(12) 



(13) 



This is an important result. However, to realize the main 
purpose of this paper, Eq. ( p^ must be made compatible 
with thermodynamics, in apparent conflict with the fact 
that it is seemingly at variance with Eq. (^), which re- 
sults in fact from the adoption of thermodynamical argu- 
ments. As a matter of fact, to establish a connection be- 
tween the dynamic and the thermodynamic approach in 
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the same spirit as that leading all the current attempts at 
unifying dynamics and thermodynamics, we must move 
from the canonical condition of to the microcanon- 
ical condition [^ . 

To prove this interesting conclusion, let us give details 
on the variational approach of Let us consider the 
functional form Fq(n) defined by 



(1-9) 



/oo 
|a;|'^n(a;)''da; 
-OO 

+a / U{x)dx. (14) 



This means that, according to the spirit of the technique 
of Lagrange multipliers, we plan to look for the entropy 
maximum under the constraint of keeping constant the 
v — th moment, second term on the r.h.s of Eq. (14), 



and the norm of the distribution n(a;), third term on the 
r.h.s of Eq. (|l|). The authors of Refs. made 
the choice ly — 2. This choice was dictated by the fact 
that in their picture the diffusing particle undergoes a 
collision process at regular time intervals, of duration r, 
making the velocity variable v fluctuate. Since x = vt, 
the constraint on the second moment of x is actually a 
constraint on the first moment of the energy distribution. 

To establish a connection between entropic and dy- 
namical arguments we are led to make a different choice. 
The space variable is made to change by the fluctuating 
velocity, which, as earlier pointed out, has only two pos- 
sible values, W and —W. During the motion in a laminar 
region, corresponding to the map of Eq. (|^), the velocity 
variable keeps the same value, either W or —W. Conse- 
quently the space x travelled is rigorously proportional 
to the time t spent in a given laminar region, and the 
transition probability 11(2:) is related to the waiting time 
distribution ^{t) by the crucial relation 



Il{x) = i;{x/W)/W. 



(15) 



This means that within the dynamic approach of this 
letter, the constraint has to be established on the first 
moment of n(a;), this being essentially a constraint on 
the mean time spent in a given laminar region. In other 
words, this physical condition implies i/ = 1. 

By means of the entropy maximization, carried out 
by differentiating Eq. (|l^) with respect to H, we get 
[A is a function of a and both Lagrange parameters are 
determined by the corresponding moment values): 



Uix) 



A 



[1 + /5(g- 1)1x1"] - 



(16) 



As argued by the authors of |T^, we note that the Levy 
diffusion process is obtained for iV oo from n(a;,iV) 
defined as the N-fold convolution product Il{x)*Il{x)*...* 
Il{x) H^]. We notice that it is not necessary to assume 
the point of view that a transition controlled by n(a:) 
takes place at any time step. In the dynamic case of 



the intermittent map of Eq. (^, we can also follow a 
different approach. We define the mean waiting time in 
a laminar region, using as a weight the distribution ip{t). 
We call this time T and we assume it to be finite. This 
is possible only if a > 1, and this is the reason why our 
dynamic approach is confined to the interval 1 < a < 2. 
Then we consider a time t » T. It is evident that for 
t —>■ oo, the distribution p{x,t) becomes very close to 
the N-fold convolution product n(a:) * n(a:) * ... * n(a;) 
with N given by the integer approximating the quantity 
t/T. As done by the authors of we apply the Levy- 
Gnedenko generalized central limit theorem p6t and we 
get a Levy diffusion with the Levy index a determined by 
the condition that the transition probability Il{x) fulfils 
the asymptotic limit 



lim n(a;) oc 



1 



\l+a ■ 



Comparing (|17|) to (l6h, we are led to 



1 



(17) 



(18) 



Note that setting ly — 1 we recover the result of the 
dynamic treatment of Eq. (p^), while the choice v — 2 
yields Eq. (^. As earlier pointed out, in the original 
work of p5[ | , the assumption was made that the energy of 
the diffusing particle fluctuates as a result of the interac- 
tion with a bath. This physical condition corresponds to 
that of Ref. ||2^, a work based on the joint use of Tsallis 
generalized thermodynamics ||l^ , and Nose-Hoover (NH) 
dynamics p^-|3l[| . In the case of ordinary thermodynam- 
ics, the NH dynamics generates, with the minimum num- 
ber of auxiliary equations, the canonical equilibrium dis- 
tribution. The work of |2^] applies the NH dynamics to 
the case of non ordinary statistical mechanics, so as to 
make energy fluctuate as in the picture adopted by Ale- 
many and Zanette |Q. 

A final issue to discuss is now the thermodynamical 
significance of the generalized version of the Kolmogorov 
entropy, namely the trajectory entropy. We note that the 
function 6{t) of Eq. (|ll|) fits very well the corresponding 
prediction of Eq. (|9). We have also seen that the func- 
tion ip{t) of Eq. (Q) agrees, through Eq. (p^, with the 
maximization of the Tsallis entropy, namely the distribu- 
tion entropy. From the comparison of Eq. (^ with Eq. 
(O) we obtain 



V>(t) = KS{~aiyo)t) 



{t > 0), 



(19) 



with cr(yo) — d^~^/yo~^ ^'^'^ = ad^^^ . A connection 
of this kind, between ip{t) and 5(t), with a independent 
of yoj is expected on an intuitive ground. In fact, the 
lifetime of a laminar region is determined by the time 
making the distance between two trajectories, initially 
very close to one another, infinitely large and so exceed- 
ing the size of the laminar region. Note however that 
'4'{t) is independent of yo and that 6{t) depends on it. 



3 



This is not surprising, since Kq^ as the conventional Kol- 
mogorov entropy Ki, is a trajectory property. To make 
Kq compatible with thermodynamics, a process forcing 
the single trajectory lose its dependence on its initial con- 
dition should be considered. The adoption of the distri- 
bution entropy seems to be more appropriate, to estab- 
lish a connection between dynamics and mechanics, as a 
consequence of the fact that it explicitly involves a sta- 
tistical average over many distinct trajectories, in accor- 
dance with the fact that also the dynamic derivation of 
the function ^{t) is obtained by making averages over 
the trajectories injected with uniform probability into a 
given laminar region. This explains why the function 
ili{t) can be successfully determined using the entropic 
argument, namely from Eq.(^|) through Eq. (|l5|), with 
no dependence on the initial conditions. 

We conclude noticing that the microcanonical treat- 
ment of this letter implies [ p5[ the mechanical energy to 
be kept fixed |^2|. The condition with the mechanical 
energy fluctuating, even if this is realized according to 
non-standard statistical mechanical prescriptions |l5[ , is 
equivalent to introducing probabilistic arguments, into a 
picture that, in principle, should remain mechanical at 
any step. In the present treatment the adoption of a 
stochastic assumption is confined to the chaotic transi- 
tion from one laminar region to the other, and it is made 
apparently unnecessary by the adoption of distributions 
rather than trajectories. 
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